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ABSTRACT 

■ Theoretical and observational investigations have indicated that the abundance of car- 

\ bon monoxide (CO) is very sensitive to intrinsic properties of the gaseous medium, 

such as density, metallicity, and the background radiation field. CO observations are 
often employed to study the properties of molecular clouds (MCs), such as mass, mor- 
phology, and kinematics. It is thus important to understand how well CO traces the 
Q ■ total mass, which in MCs is predominantly due to molecular hydrogen (H2). Recent 

Jjj ' hydrodynamic simulations by Glover & Mac Low have explicitly followed the for- 

mation and destruction of molecules in model MCs under varying conditions. These 
models have confirmed that CO formation strongly depends on the cloud properties. 
Conversely, the formation of H2 is primarily determined by the amount of time avail- 
. able for its formation. We apply radiative transfer calculations to these MC models 

J> 1 in order to investigate the properties of CO line emission. We focus on integrated 

ON ! CO (J=l-0) intensities emerging from individual clouds, including its relationship to 

the total, H2, and CO column densities, as well as the "X factor," the ratio of H2 
column density to CO intensity. Models with high CO abundances have a threshold 
CO intensity «65 K km s _1 at sufficiently large extinctions (or column densities). 
Clouds with low CO abundances show no such intensity thresholds. The distribution 
of total and H2 column densities are well described as log-normal functions, though 
the distributions of CO intensities and column densities are usually not log-normal. In 
general, the probability distribution functions of the integrated intensity do not follow 
the distribution functions of CO column densities. In the model with Milky Way-like 
conditions, the X factor is in agreement with the near constant value determined from 
' observations. In clouds with lower metallicity, lower density, or a higher background 

H \ UV radiation field, the CO abundances are in general lower, and hence the X factor 

can vary appreciably - sometimes by up to 4 orders of magnitude. In models with high 
densities, the CO line is fully saturated, so that the X factor is directly proportional 
to the molecular column density. 

Key words: ISM: clouds - ISM: molecules - ISM: structure - methods: numerical - 
stars: formation 



1 INTRODUCTION 

Since stars form almost exclusively in clouds which are pre- 
dominantly molecular, much effort in current star formation 
research is focused on understanding the physics and chem- 
istry of molecular clouds (MCs) . Though MCs are primarily 
composed of molecular hydrogen (H2), due to the lack of a 
dipole moment and the unsuitable conditions within MCs 
to excite its rotational transitions, H2 is difficult to observe 
directly. Helium, accounting for ~ 10% of the cloud mass, is 
also very difficult to detect. Carbon monoxide (CO), which 
is the second most abundant molecular species in MCs, has 



a dipole moment with rotational transitions that are easily 
excited at typical MC temperatures (10 - 100 K) and densi- 
ties ( 100 cm -3 ). CO observations are thus often employed 
to investigate MC properties. 

The observed CO intensity, I (CO), which is often 
expressed as a integrated "brightness temperature" Wco 
(hereafter W), is considered to be a good tracer of the col- 
umn density of molecular hydrogen JVh 2 • Namely, Nh 2 is 
estimated from CO obser vations through a constant U X fac- 
tor" fe.g. |Pickmanlll97Sl) : 
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x = 



Nh 2 

w 



1 KT 1 km" 



(1) 



CO observations of Galactic clouds have resulted in es- 
timates of X w few xlO 20 cm" 2 K" 1 km -1 s (here- 
after X Ga i, e.g.lSolomon et alJll987l : lYoung fc Scovillelll99ll : 
iDame et al.ll2Q0lh . Observations of diffuse gas in the Galaxy 
have also resulted in similar estim ates of the X factor (e.g. 
iPolk et al.|[i"98Sl ; iLiszt et alj|20ich . 

However, extragalactic observations of systems with dif- 
ferent physical characteristics, such as metallicity or back- 
ground UV radiation, have found variations in the X factor. 
Interestingly, observational investigations employing differ- 
ent methodologies have resulted in vastly discrepant esti- 
mates of the X factor. Fo r example, for the n earby Small 
Magellanic Cloud (SMC), iBolatto et ail IMPOST } measure a 
value ~Xg & \, assuming virialized clouds and using the CO 
linewidths to estimate cloud masses, and thereby Nn 2 • In- 
dependent dust and gas based observations, on the other 
hand, suggest extended regions containing molecular mate- 
rial with little or no CO emission, resulting in an X factor of 
up to ~100 X Ga i jlsraell 19971 ; iRubio et al.ll2004l ; iLerov et all 
I2007l.l2009h . 

The 12 CO (J=l-0fl line is optically thick in most molec- 
ular clouds, and so CO (J=l-0) observations are known not 
to provide direct information about the total CO mass or 
column density -/Vco- Observations have shown a satura- 
tion of CO intensities at sufficiently high extinctions (e.g. 
lLombardi et all 120061 ; iPineda et all 120081 ). Consequently, 
only CO intensities below the saturation threshold are con- 
sidered in evaluations of the X factor. One of the goals of 
this work is therefore to assess how well CO observations 
can trace the true CO distribution within molecular clouds. 

Since the formation of CO is sensitive to the amount 
of carbon and oxygen available in the ISM, CO abun- 
dance s are expected to be lower in low er metallicity sys- 
tems l|Malonev fc Black|[T988l ; llsraellll99^ . Additionally, the 
strength of the background ultraviolet (UV) radiation field, 
responsible for photodissociating mol ecules, also plays a role 
in re gulating the CO abundances (|van Dishoeck fc Blackl 
1988). These processes should in turn lead to variations 
in the X factor, with a larger value in metal poor sys- 
tems and/or where the UV radiation is higher than in 
the Milky Way. Indeed, when independent measures of MC 
masses are combined with CO observations, the X factor has 
been found to be larger in lo w metallicity external galaxies , 
such as the LMC and SMC l|lsrael et al.lll986l ; Ilsraellll997l ; 
ILerov et al]|2009D . 

As turbulence is now considered an important aspect 
of the dynamics in molecular clouds and star formation 
jMac Low fc Klessenl 120041 ; iBallesteros-Paredes et al.ll2007l ; 
iMcKee fc Ostrikeil I2007T and references therein), its role 
in influencing mole c ule form ation need s to be understood. 



Glover fc Mac Low! ( |2007bh (see also iGlover fc Mac Low! 
2007al ), showed that the formation timescale of H2 is only 



a few Myr in turbulent MC models with densities compara- 
ble to those in observed clouds. Expanding on this work, 
through implementation of more extensive chemistry, we 
can now follow the fo rmation of CO in models of turbu- 
lent molecular clouds (|Glover et alj I2OI0I . hereafter Paper 



1 We will hereafter refer to 12 CO simply as CO. 



I). Subsequentlv. [G~lover fc Mac Low! |2010T l. hereafter Paper 
II, analysed the global properties of the X factor in different 
MC models, using spatially averaged quantities. They found 
that the mean extinction, or total gas column density, pri- 
marily determines the X factor. 

Here, we explore the observational consequences of the 
variation in CO abundances which arise due to differences in 
MC properties described in Papers I-II. Our main goal is to 
understand the impact of CO abundance variations within 
individual MCs on the emerging integrated CO (J=l-0) in- 
tensity. We consider a suite of models with different condi- 
tions, namely metallicity, density, and background UV radi- 
ation field, representing various environments. In the next 
section we discuss our modeling method, including a brief 
overview of the magnetohydrodynamic models with chem- 
istry, and a description of the radiative transfer calculations. 
In Section [3] we present our results of the comparison of CO 
intensities with H2 and CO column densities, and charac- 
teristics of the the X factor. We discuss our results and 
compare them to observational investigations in Section [4] 
We conclude with a summary in Section [5] 



2 MODELING METHOD 

To carry out our investigation of CO emission from molec- 
ular clouds, we apply line radiative transfer calculations to 
magneto-hydrodynamic (MHD) models of MCs. In this sec- 
tion, we give a brief description of the MHD models, and 
focus on the radiative transfer calculations. We refer the 
reader to Papers I and II for a more extensive description 
of the MHD and chemical modeling method, as well as the 
analysis of various simulation runs. 



2.1 Modeling Molecular Clouds: MHD and 
Chemistry 

The simulations of the model MCs track the evolution of 
gas using a modified version of the Z EUS-MP MHD code 
|Normanl200(il ; IStone fc Normanll992al lbh. Gas with initially 
uniform density in a 20 pc box with periodic boundary con- 
ditions is driven with a turbulent velocity field, with uniform 
power between wavenumbers 1 ^ k sC 2. Gas self-gravity, 
which would cause sufficiently overdense regions to collapse, 
is not considered in the calculation. The simulation includes 
magnetic fields, with initial field lines oriented parallel to 
the 5-axis with a field strength of 1.95 /J.G. The gas has an 
initial temperature of 60 K, but quickly settles to thermal 
equilibrium. 

To track the chemical evolution of the gas, which has 
constant metallicity, a treatment of hydrogen, oxygen, and 
carbon chemistry is included in the numerical algorithm. 
This consists of 218 reactions of 32 chemical species, which 
are coupled with the thermodynamics. A background UV 
radiation field, which is responsible for photodissociation, is 
treated through th e six-ra y approximation as described in 
IGlover fc Mac Low! (|2007al ). In our investigation of emission 
from CO molecules, we consider one snapshot of these sim- 
ulations at a late time > 5 Myr, after which the simulation 
has reached a statistically steady state. 
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2.1.1 GMC Model Parameters 

The numerical modeling described above follows the com- 
bined effects of turbulence, magnetic fields, thermodynam- 
ics, and chemical evolution as structures such as filaments 
and dense cores typically found in MCs forms out of the gas 
(see Paper I for more details). In our investigation of CO 
emission emerging from the MC simulations, we consider a 
suit of models designed to represent various astrophysical 
environments. The relevant parameters of each model are 
listed in Table [T] Column 1 shows the name of each run. 
The main user defined parameters are the initial density no, 
metallicity Z, and background UV radiation field Go, in- 
dicated in Columns 2-4, respectively. We assume that the 
dust-to-gas ratio is directly proportional to the gas-phase 
metallicity and do not vary these quantities independently. 
Column 5 shows the numerical resolution; as noted in Paper 
II, simulations with identical initial conditions but resolu- 
tions of 256 3 and 128 3 produce very similar results0 The last 
column lists the representative environment corresponding 
to the simulated cloud: a high density cloud found in the 
galactic cente r (nlOOO), a typical Milky Way cloud (n300, 
lFerriere|[200ll ). clouds in low metallicity systems like the 
LMC or SMC (n300-Z03 and n300-Z01), a low density cloud 
in a dwarf galaxy (nlOO), and clouds in weak and strong 
starbursts (n300-UV10 and n300-UV100, respectively). 

2.2 Radiative Transfer Method 

2.2.1 Radiative Transfer Overview 

We use the radiative transfer code RADMC-3D (Dullemond 
et al. in preparation) to model the molecular line emission of 
the MHD MC models. RADMC-3D is a 3-dimensional code 
that performs dust and/or line radiative transfer on Carte- 
sian or spherical grids (including adaptive mesh refinement) . 
In this work, since we are only interested in CO molecular 
line emission along a chosen direction, we use the ray tracing 
capability of RADMC-3D. 

One of the primary challenges in line radiative trans- 
fer is to solve for the population levels of the molecular (or 
atomic) species under consideration. The occupation of a 
given (rotational/vibrational) energy level of a molecule de- 
pends on the incident radiation field, as well as the collision 
properties (e.g. frequency) with other atoms or molecules, 
both of which act as excitation or de-excitation mechanisms. 
In statistical equilibrium, the relative population of level i, 
fi, is governed by the equation of detailed balance: 



y ] [.fjAji + (fjBji fiBij)Jji] — 

j<i 



(2) 



collisional rate for a transition from level i to level j. The 

2 We also find few differences in the results from the radiative 
transfer calculations applied on those models with different reso- 
lutions but otherwise identical initial conditions. 



collisional rate is dependent on the rate coefficient K%j and 
the density of the collisional partner n co i: Cy = n co iKij. In 
MCs, the main collisional partner of CO is H2. We use the 
rate coefficients for collisional excitation and de-excitations 
of CO by H2 tabulated and fr eely available from the Leiden 
databas e (ISchoier et alj|20 05). These are based on calcula- 
tions bv lYang et al. I (|201C| yr We neglect the effect of collisions 
with partners other than H2. We justify this by noting that 
almost all of the CO in our simulations is found in regions 
of the gas that are dominated by H2, rather than atomic 
hydrogen. 

The influence of radiation on setting the level popula- 
tions is captured by the first two summations in Equation 
[2] which include the mean integrated intensity Jij of the ra- 
diation field in the line corresponding to the transition from 
i to j. The constants are Aij, the Einstein coefficient for 
spontaneous emission for a transition from level i to level 
j, Bij, the Einstein coefficient for stimulated emission from 
i to level j, and Bji, the corresponding coefficient for ab- 
sorption (also given by the Leiden database). Equation [2] is 
coupled with the equation of radiative transfer 

^ = -/„+<?„, (3) 

where I v is the specific intensity, S v the source function, and 
t v is the optical depth. The coupling between Equations [2] 
and [3] occurs through the dependence of the source function 
[S v = Sij) on the relative population levels: 



fjBji fiBij 
as well as the dependence of the J on I: 



Jij — ^ j lij(f)ijdVl : 



(4) 



(5) 



where the integral is taken over all solid angles Q. The nor- 
malized profile function (j>ij determines the emission and ab- 
sorbtion probability of a photon with frequency u, due to a 
line with rest frequency Vij. For a photon propagating in a 
direction n through a medium with velocity v, 



4>ij{v) = 



■. exp 



c 2 (u — Vij — v • huij/cy 



(6) 



where c is the light speed. The thermal and microturbulent 
broadening of the line is accounted for through 

a 2 = Vmtrb + 2fcsT/m mo i, (7) 

where fcs is the Boltzmann constant, m mo ; is the molecu- 
lar (or in the case of atomic lines, atomic) mass, T is the 
temperature, and v m trb is the microturbulent velocity of the 
medium. 

As Equations [2] - [5] indicate, the amount of radiation 
emitted from molecules in a given location is dependent 
on the level populations, which themselves depend on the 
amount of incident radiation at that location. Solving this 
problem numerically can be computationally expensive. In 
certain situations, however, suitable approximations can be 
made which considerably reduce the computational costs. 

One approximation is that of local thermodynamic equi- 
librium (LTE). In gaseous systems, when collisional pro- 
cesses dominate line emission, the temperature is the only 
parameter that is required to calculate the population lev- 
els, through the use of a partition function. Thus, in regions 
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Table 1. List of Simulations 



ID 


no (cm 3 ) 


Z 


Go (2.7xl0~ 3 org cm" 2 s" 1 ) 


Resolution 


Representative Environment 


nlOOO 


1000 


1.0 


1.0 


128 3 


high-density cloud (galactic center) 


n300 


300 


1.0 


1.0 


256 3 


Milky Way cloud 


n300-Z03 


300 


0.3 


1.0 


256 3 


LMC/SMC cloud 


n300-Z01 


300 


0.1 


1.0 


256 3 


LMC/SMC cloud 


nlOO 


100 


1.0 


1.0 


256 3 


dwarf galaxy cloud 


n300-UV10 


300 


1.0 


10.0 


128 3 


weak starburst 


n300-UV100 


300 


1.0 


100.0 


128 3 


strong starburst 



with high Hb densities, employing the LTE assumption to 
model the CO (J=l-0) line is reasonable. However, in the 
MC environments considered here, much of the volume does 
not contain high densities, so the LTE approximation may 
not be suitablejj 



2.2.2 The Sobolev approximation 

For this study, we hav e implemented the Sobolev approx- 
imation l|Sobolevl Il957l : also known as the Large Velocity 
Gradient, or LVG, method) into RADMC-3D to solve for 
the population levels of an atomic or molecular species. This 
method takes advantage of large spatial variations in veloc- 
ity, as are present in turbulent molecular clouds, to define 
line escape probabilities. Effectively, the LVG method pro- 
vides a solution to the equation of detailed balance (Eqn. [2]) 
determined solely by local quantities. 

Consider a photon emitted due to a transition from 
level i to level j. Beyond a certain distance from the pho- 
ton's position of emission, due to the velocity gradient in 
the medium, the Doppler-shifted frequency associated with 
the i to j transition is sufficiently different from that of the 
incident photon. The photon cannot interact with matter at 
this position, or any subsequent position, and thus propa- 
gates freely out of the cloud. The LVG escape probability of 
the photon can be determined from the optical depth 



T dr' = 



(8) 



Other functional forms of the escape probability can be de- 
fined fo r different scenarios (e.g . slab or spherical symme- 
try, see Ivan der Tak et al.l [20071 ). The optical depth is de- 
termined by the population levels, dens ities, and velocity 
gradients dv/dr ijvan der Tak et aill2007l ): 



AijTl 



8tti/ 3 . 1.065 < \dv/dr\ > 



9i 



(9) 



where gt is the statistical weight of level i, and n is the 
total density of the molecular species (fi is the relative pop- 
ulation, so that ^2 fi = 1). We use the mean (absolute) 
velocity gradient < \dv/dr\ > across the six faces of each 



3 We have also performed radiative transfer calculations assum- 
ing LTE to obtain CO (J=l-0) line intensities for all the GMC 
models considered in this work. As expected, the LTE intensities 
only differ by a factor of a few from the LVG intensities. Had we 
only considered LTE, our overall (qualitative) conclusions would 
be unchanged. The precise values of W and the X factor, however, 
would be appreciably different. 



grid zone in the simulation. The near unity factor in the 
denominator is a correction for integration over a Gaussian 
line profile Q 

Once the optical depth is computed, the escape prob- 
ability is given by Equation [8] The local radiation field is 
then determined by 



jy = Sy(l-/3) 



(10) 



Numerically, given these expressions for r, /3, J, and the 
local velocity gradient, one can iteratively solve for the pop- 
ulation levels in Equati on [2] for more detai led description s 
of the LVG method see iMihalasI \ 19781 ) and lElitzurl j 19921 ) . 

Though the LVG method was originally formulated and 
developed for spectral line studies i n stellar atmospheres 
which have smooth velocity gradients l|Sobolevlll957l ; ICastorl 
1970; Mihalas 1978), it can also be applied for a ny environ- 
ments with significant velocity gradients. In fact. lOssen kopl 
l|l997l) showed that the Sobolov method is a very good ap- 
proximation for computing level populations in MCs, where 
the velocity gradients are generally not smooth0 



2.3 Analysis of CO Emission 

In this work, we focus on how well CO emission can trace the 
intrinsic CO column density, Nco, the H2 column density 
Nh 2 , as well as the total column density of hydrogen nuclei, 
iVtot • The intrinsic column densities can be easily calculated 
directly from the MHD simulation by integrating the CO, H2 
and hydrogen nuclei volume densities along a given axis. Ac- 
cordingly, for the radiative transfer calculation we orient the 
simulation cube such that the CO line is "observed" along 
the same axis for which the column densities are computed 
(e.g. along the f-axis). We have verified that our results are 
not sensitive to the choice of orientation. 

Besides the viewing geometry, the only other user de- 
fined parameter required for the radiation transfer calcula- 
tions is the microturbulent velocity v mtr b [see Eqns. I6l7j . 
Extrapolating the observed linewidth-size relationship (e.g. 



4 We have verified that our LVG implementation is accurate by 
comparing the derived population levels of various simple mod- 
els against those provide d by the RADEX online LVG caculator 
[Ivan der Tak et alj|2007fl . and that the LVG solution at high den- 
sities matches the LTE solution. 

5 A photon determined to escape through an LVG calculation 
may be absorbed in a distant region for a system with a stochastic 
velocity distribution. We do not consider such events, as they 
are not expected to occur frequently in systems that are highly 
supersonic. 
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iLarsonl Il98lh down to the resolution of the MHD simula- 
tions ~ 0.1 pc, appropriate microturbulent velocities are in 
the range ~ 0.2 - 0.7 km s _1 . We have explored this range 
in v m trb> and have found that the results are insensitive to 
the particular choice of v m t r b- In the analysis presented here, 
Vmtrb is set t° 0.5 km s _1 . 

In our discussion, we will sometimes express ATtot as 
an extinction Av, in order to allow for direct comparison 
with observational analyses. Observers are often required to 
employ indirect measures of A^ot, since H2, a major con- 
stituent of TVtot in MCs, is difficult to observe directly. Ex- 
tinction measurements provide estimates of the total amount 
of dust along the line of sight. Using a "reddening law," 
(and an assumption for the dust-to-gas ratio), the total 
gaseous c olumn follows direc tly from the amount of extinc- 
tion Av (Boh lin et al.lfl97Sl) . In many Galactic molecular 
clouds, nearly all the hydrogen is molecular, so Nn 2 is di- 
rectly proportional to Ay. In other environments, however, 
such as those with lower metallicity or higher UV fields (e.g. 
Models n300-Z01 and n300-UV2), there may be significant 
amounts of atomic hydrogen, so Ntot, and hence the extinc- 
tion, is dependent on both the molecular and atomic column 
densities. To allow for straightforward comparison between 
models and observations, we use a simple conversion be- 
tween iVtot and Ay: 

Av= 1.87xK£cm-a (^) ' 

where Z is the metallicity of the gas. Thus, in the com- 
parison of CO emission with intrinsic cloud properties, any 
discussion involving Av can be directly translated into total 
column density. 

The radiative transfer calculations produce spectral 
(position-position- velocity, or PPV) cubes of the CO (J=l- 
0) line. The 3D cube indicates the intensity I v in a given fre- 
quency or velocity channel of width dv at each 2D position. 
We choose a sufficiently large range in frequencies such that 
all (line-of-sight) velocities in the simulation are detected, 
so that all emission from the model MC is "observed." 

In the comparison of CO intensities with intrinsic col- 
umn densities, as well as the analysis of the X factor, the 
quantity of interest is the velocity integrated intensity, which 
is simply the PPV cube integrated over its velocity axis. 
The intensity /„, which has units of erg s _1 cm -2 Hz -1 
ster -1 , can be expressed as the Planck function evaluated 
at a "brightness temperature" Tb, B v (Tb)- Since the CO 
(J=l-0) line is located in the Rayleigh- Jeans part of the 
spectrum, l v oc Tb- The intensity of CO line emission is 
thus often conveyed in Tb units. We follow this convention 
and express the intensity as a velocity integrated brightness 
temperature: 

W = -^-{c/vf / I v dv (Kkms" 1 ). (12) 
2k B J 

This integrated intensity is thus a measure of total CO emis- 
sion along the line of sight. W is computed at all positions 
yielding a 2D map, which can be used with the 2D map of 
JVh 2 to obtain the X factor through Equation [T] 



3 RESULTS 

3.1 CO Intensities and Column Densities 

To begin our investigation of CO emission, we assess the 
correlation of velocity integrated CO intensities W with the 
total column density Ntot of the model MCs. Figure [1] shows 
W as a function of the N to t for four of the simulations listed 
in Table Q] The bottom abscissa shows Ntot , while the top 
abscissa shows the corresponding Av (see Eqn. Ill[) . 

A number of features are readily apparent in Figure [T] 
First, there is a general trend of increasing intensity with 
increasing column density, though the slopes differ between 
the various simulations. For the high-density run (nlOOO) 
shown in Figure [T};, the CO intensities do not demonstrate 
any clear trend with increasing Av- Rather, the vast ma- 
jority of the intensities reach a threshold value of ~65 K 
km s _1 . This saturation of CO intensities is expected to oc- 
cur at high densities since the CO line becomes optically 
thick. Saturation is also found at the highest extinctions in 
the Milky Way simulation (n300, FigQJ,). Indeed, satura- 
tion at high CO intensit ies has been observed in MCs in the 
solar neighborhood (e.g. lLombardi et al. l l2006l ; |Pineda~ et al.l 
2008). Though model n300 qualitatively reproduces the ob- 
served trends of increasing CO intensity with increasing den- 
sity, up to a threshold value, a detailed quantitive compari- 
son is not appropriate here; the precise slope, scatter, min- 
imum, and threshold intensity are dependent on additional 
physics not included in our models, such as additional heat- 
ing due to stars, outflows, and supernovae. 

Figure [Tp and d show that for clouds with lower metal- 
licities (Model n300-Z03) or densities (Model nlOO), there 
is a wider distribution of CO intensities at low extinctions. 
The saturation of the CO line is not easily evident compared 
to Model nlOOO. Since CO is easily photo-dissociated in low 
metallicity or low density systems due to insufficient (dust 
and self-) shielding, the CO column densities do not reach 
high enough values for the saturation of the CO line to occur. 
Observations have indeed shown that the CO intensities are 
generally lower in extragalactic systems with lower metallic- 
ities than that of the Milky Way. Further, there is no evi- 
dence of saturation in the CO line from di ffuse sources such 
as the Polaris Flare, LMC and SMC (e.g. iHeithausen et all 
1 19931 ; Ilsraellll997l ; iLerov et alj|2009h . 

Figure [2] shows the probability distribution functions 
(PDFs) of JVtot (or A v ), Nn 2 , N C o, and W for five mod- 
els. The Ntot PDFs are offset to larger values compared with 
the A?h 2 PDFs. The offset for the Milky Way (n300) and the 
high density cloud (nlOOO) can be wholly attributed to the 
difference in Ntot, which is a measure of the total number 
of hydrogen nuclei, and A% 2 , a measure of the number of 
molecular hydrogen nuclei. For the other models, there re- 
mains an offset even after accounting for the factor of 2 dif- 
ference between Ntot and Nh 2 , because significant amounts 
hydrogen remains in atomic form (see Table 2 of Paper II) . 

The Ntot and A% 2 PDFs of all simulations can be well 
described as log-normal functions. Log-normal column den- 
sity PDFs should b e expected from these compressible tur- 
bulent simulations (IVazquez-Semadenilll994l; iKle sscn 2000; 
lOstriker etail bOOll - lTi et alj|2003l : iFederrath et alj|20ich . 
The simulation with conditions similar to the Milky Way 
(n300) and with the highest density (nlOOO) produce Nco 
and W PDFs which come closest to being lognormal. In 
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general, however, the CO column densities and intensities 
do not show a log-normal PDF. We defer the discussion of 
the best fit Gaussians shown for these two models to Section 

ED 

Another clear trend in Figure [2] is that the shapes of 
the W PDFs can differ significantly from the PDFs of un- 
derlying Aco- The extent of the scales of W, shown on the 
top abscissas of the right hand plots in Figure [2] is equiv- 
alent to the extent of the scales in Aco, which appear on 
the bottom abscissas. At low densities where gas is optically 
thin, CO emission should trace the CO density. Thus, the 
absolute scales are chosen such that low W and Aco values 
are matched as best as possible. For models n300 (Figs[2Ji-b) 
and nlOOO (Figs [2]>f) , the W PDFs have steep increases, 
resulting in peaks that do not correspond to the peaks in the 
A CO PDFs. The PDFs of the other simulations, n300-Z03, 
n300-Z01, and nlOO, show comparable gradients at low W 
and Aco , but at high values the W PDFs drop more rapidly 
than Aco- 

In those PDFs showing a steep decline in W, at high 
CO column densities there is little or no corresponding CO 
emission, even though there is good correspondence at lower 
values (due in part to the choice of the W and Aco axis 



ranges). PDFs showing a much steeper gradient beyond the 
W peak compared to the Aco peak can be attributed to sat- 
uration of the CO line: as a result of saturation, regions with 
very high A^co may have similar intensities to regions with 
lower (but still high) Aco, producing a PDF which appears 
to have a "piled-up" profile at some high W. Such PDF 
shapes are clearly evident in the moderately low metallicity 
(n300-Z03) and very low density (nlOO) cloud simulations. 
For the simulation with very low metallicity (n300-Z01) for 
which the W and Aco PDFs appear to be correlated, there 
is no evidence of a saturation feature. 

We have shown that there is no simple correlation be- 
tween W and Aco, focusing on models with different metal- 
licities and densities. Of course, the background radiation 
field also affects the formation of CO, due to UV photodis- 
sociation |van Dishoeck fc Black! [l98Sl , Paper II). Figure [3] 
shows the Ah 2 , Aco, and W histograms from models n300, 
n300-UV10, and n300-UV100 (see Table D). These models 
only differ in their background UV radiation field strengths, 
with values of Go = 1, 10, and 100 times the Galactic es- 
timate (2.7xl0 -3 erg cm -2 s _1 ). All simulations have sim- 
ilar log-normal Nh 2 distributions. Simulations with higher 
UV fields have more extended iVco PDFs, and thereby low- 
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ered maxima in the Nco probabilities. This occurs because 
stronger UV radiation is able to penetrate further into the 
cloud to dissociate CO, resulting in more regions with low 
CO densities. Similarly, the W distributions are more ex- 
tended for models with brighter UV backgrounds, with fewer 
regions with high W values. These results simply follow the 
trend previously discussed: models with a wide range in CO 
abundances demonstrate a lack of correlation between W 
and iVco PDFs. The lack of correlation in the PDFs primar- 
ily occur at the highest densities and intensities (with the 
exception of n300-Z01, but see last paragraph of this subsec- 
tion), though at lower densities gas is optically thin result- 
ing in a correlation between W and Nco ■ Clouds with lower 
metallicity, lower overall density, higher UV background, or 
a combination of these factors have much more distributed 
CO abundances, resulting in greater discrepancies between 
the W and iVco PDFs. Keeping this general trend in mind, 
in our subsequent analysis we will focus on the models with 
Go = 1. 

We now turn our attention to the W and Nco proper- 
ties of simulation n300-Z03 (Figs.[2j;-d). There is a relatively 
large spread in Nco, which also has a clear double peak, at 
log(JVco) ~ 14.5 and 16.5 (with similar probabilities). Fur- 
ther, though the W PDF appears to correlate well with the 
Nco PDF at low CO density (or low intensity) values, the 
peak in W does not correlate well with the peak in Nco- A 
closer inspection of the Nco and W images for Model n300- 
Z03, shown in Figure [4ji-b, further reveals the correlation 
at low densities, and lack of correlation at high densities. In 
the low intensity range —3 ^ log(W) ^ — 1, marked by solid 
lines in Figure [2}i, an increase in W corresponds to a similar 
increase in iVco in the 12 ^ log(A?co) ^ 14 range. These 
W and Nco values are indicated by solid contours in Figure 
[4] showing excellent correspondence. However, at higher in- 
tensities, W does not faithfully trace the high iVco regions. 
The PDF of intensities at log(W) 1.5, marked by a dashed 
line in Figure [2}i, decreases sharply, whereas log(iVco) con- 
tinues to increase. The dashed contours in Figure [3] mark 
the corresponding regions in the iVco and W maps. In some 
of the filaments (e.g. near the bottom), the contours mark 
similar regions. However, in the most dense regions (e.g. the 
large over-dense region in the top right), the area within 
the W contour is noticeably smaller than that encompassed 
in the iVco map. The detailed morphology in the region 
encompassed by the dashed contour also shows stark differ- 
ences between the observed CO intensities and intrinsic CO 
column densities. These differences at high densities are a 
direct consequence of the optically thick nature of the CO 
line. 

Similar discrepancies are present in a Nco - W compar- 
ison of all the models. In some cases, such as model cloud 
with very low metallicity, n300-Z01 shown in Figure [2], the 
apparent correlation in PDFs at low densities does not trans- 
late into a correlation in the 2D maps. As stated, the Nco 
and W PDFs were plotted with abscissas chosen to best 
match the PDFs at low densities and intensities, and so any 
apparent correlation may only be coincidental. 

In summary, we have shown that at low densities, CO 
emission traces the distribution of CO molecules well, since 
the CO line is optically thin. However, at high densities, the 
line becomes optically thick, resulting in an intensity thresh- 
old. Thus, CO intensity does not neatly trace the highest 



density gas. Further, there is generally no simple correla- 
tion in the distribution of CO and H2 molecules, since H2 
can better shelf-shield against photodissociation. Taken to- 
gether, the lack of any general correspondence between W 
and iVco, as well as iVco and A r H 2 , will certainly affect fur- 
ther analysis of CO observations, such as attempts at relat- 
ing W to iV H2 • 

3.2 The X Factor 

In Paper II, the X factor (Eq.[TJ from the models considered 
in this work was found to vary at low mean extinctions Av 
< 3. At higher mean extinctions, the X factor remained 
constant. In that work, the CO intensities were estimated 
by assuming local thermodynamic equilibrium (LTE), and 
by averaging relevant quantities over whole clouds. In the 
previous section, we demonstrated that CO observations do 
not neatly trace CO column density. A logical expectation 
from this result is that observed CO intensities would not be 
a good tracer of H2 number density, leading to a variation 
in the X factor as found in Paper II. In this section, we go 
beyond the analysis of Paper II to investigate the X factor 
using the results from the 3D radiative transfer calculations 
to estimate CO line intensities. Here, we are also able to 
assess variations in the X factor on a position-by-position 
basis in maps of the individual clouds. 

Figure [5] shows the X factor from each position in the 
2D integrated intensity images of simulations n300, n300- 
Z03, nlOOO, and nlOO, along with the PDF of the X factor. 
Each point indicates the relationship between iVco and Nn 2 , 
and the color identifies its value of the X factor, calculated 
from Equation [1] 

For the Milky Way MC model (n300, Fig.[5Ji), in >90% 
of positions the X factor is within the range 9x 10 19 — 5x 10 20 
cm~ 2 K _1 km -1 s. Figures [2^,-b show that the quantities 
Nh 2 , Nco, and W all have similar (Gaussian- like) shapes, 
which ultimately result in a limited range in the distri- 
bution of the X factor. The range in the derived X fac- 
tor for this model is in agreement wit h estimates from ob- 
servations of Milky Way clouds (e . g ISolomon et al.l 1 19871 : 



lYoung fc Scovillelll99ll : iDame et all 



20011 ). This agreement 



gives us confidence that the simulations, including the chem- 
ical networks as well as the radiative transfer calculations, 
are sufficiently reliable and can be utilized to probe condi- 
tions different from the Milky Way. 

The distribution of the X factor increases substantially 
when the observed CO emission is poorly correlated with 
the CO column densities, as in Model n300-Z03 (Fig. [Sb). 
In this model, the X factor distribution spans a wide range: 
only < 25% have values in the range 9x 10 19 - 5x 10 20 , and 
>30% have values > 10 22 cm -2 K" 1 km -1 s. 

Figure Hfc and d shows Nb_ 2 and the X factor for each 
line-of-sight of Model n300-Z03. In the highest density re- 
gions such as the large complex at the top right-hand side 
of CO or H2 column density maps (Figs. [4^ and c) , the X 
factor has a limited range, within one order of magnitude 
(10 20 — 10 21 cm" 2 K _1 km -1 s). In those regions, though 
Nco can vary by 2 or more orders of magnitude, W only 
varies by « 1 order of magnitude. Since the X factor di- 
rectly depends on W and only indirectly on iVco, the X 
factor also only falls into a limited range. 

Positions with the largest X factors correspond to the 
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lowest Nh 2 regions, as well as low A^co and W regions. These 
are the regions where CO is most affected by photodissocia- 
tion. Since the amount of photodissociation depends on the 
the "effective" column density in each location of the 3D 
simulation volume, regions with similar JVh 2 can have very 
different ATco values (see also Paper I and II), as evident in 
Figure [5]d: at low to intermediate H2 densities 10 21 -/Vh 2 

10 22 cm -2 , there is a wide range of iVco for a given 
JVh 2 • Since the X factor (indirectly) depends on A^co , at 
such densities the X factor also takes on a wide range. For 
instance, at Nu 2 — 5 x 10 21 cm -2 , the X factor varies from 
~ 10 20 to 10 23 cm -2 K _1 km -1 s. Evidently, the X factor 
can have a wide distribution within a MC, even for regions 
with identical molecular column densities. This is a conse- 
quence of the combination of a large distribution of iVco for 
a given Nn 2 , as well as the lack of simple correlation between 
W and JVco due to the optically thick nature of CO. 

Figure [5]i shows that there can be a wide range in the 
X factor even in very low density regions. For this model (no 
= 100 cm" 2 and Z = Z Q ), much of the gas has Nh 2 ~ 10 21 
cm -2 . Unlike Model n300 (in Fig. [5^), there is a very wide 
distribution in the X factor in the range 10 20 < JVh 2 £ 10 21 
cm -2 . This model also differs from Model n300-Z03 (in Fig. 
[5p), showing a much larger distribution in Nco and the X 
factor for a given Ah 2 at Ah 2 ~ 2 x 10 21 cm" 2 . The multiple 
peaks in the X factor distributions can be understood from 
the W and Nco distributions shown in Figure [2ji. Model 
nlOO has a large distribution of Nco', the W distribution is 



similar, with multiple peaks (which are not aligned with the 
Nco peaks), contributing to the wide range and multiple 
peaks of the X factor shown in Figure [SJl. 

At very high densities, more self shielding allows for 
CO formation to be more effective. As Figure [Sfc shows for 
Model nlOOO, Nco is a much better tracer of ATh 2 • However, 
though A'h 2 for this model reaches values 3 times greater 
than Model n300 (Fig.Hk), the distribution in the X factors 
of both models is very similar. In Model nlOOO, the X factor 
increases with increasing (CO or H2) density, a trend oppo- 
site to that apparent for the other models shown in Figure 

El 

The correlation between the X factor and Ah 2 is shown 
in Figure El Model n300 shows that at A% 2 £ 10 21 cm -2 , 
the X factor is ~ constant (with a mean value 2 x 10 20 
cm" 2 K _1 km" 1 s). Only at the highest densities (log(A r H 2 ) 
£ 22.0) do the effects of saturation become clearly evident. 
On the other hand, Model nlOOO (shown in Fig. [BJ;) shows 
an increase in the X factor at all ATh 2 . The line shows the 
X = Ah 2 /(70 K km s" 1 ) relationship, which is simply the 
H2 column density divided by the mean CO intensity from 
this model. This line is a very good fit at log(A r H 2 ) ~ 22, 
indicating that at most densities the CO line emission is 
completely saturated. 

The models with lower metallicity and density (Figs.[6p 
and [SJl) show a decreasing X factor with increasing density 
(as can also be seen in Fig. [5]). In this regime, an increase in 
the molecular density is associated with a greater increase 
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a) log(N C0 ) (cm -2 ) b) log(W) (K km s" 1 ) 
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Figure 4. Images of a) Nqo, b) W, c) Nn 2 , and d) the X factor of Model n300-Z03. Each side has a length of 20 pc. In a) and b) solid 
contours indicate log(A r co) = 12, 14 and log(W / )=-3, -1; dashed contours are log(7Vco)=16.5 and log(W)=1.5 (see text and Fig.[2ji). 



in CO line intensity, resulting in a decrease in the X factor. 
Such a trend is also present at the lowest densities of Model 
n300 (Fig. Hi). 

Given that W can have rather different distributions 
than -/Vco, it comes as no surprise that the X factor does 
not neatly follow from iVco ■ For clouds with high CO abun- 
dances, such as models n300 and nlOOO, Nco is well corre- 
lated with Nb_ 2 , resulting in a limited range in the distribu- 
tion of the X factor. We have seen that CO line saturation 
may not be discernible solely from the X factor distribu- 
tions, as evident from the similarity of the X factor (nor- 
mal) distributions of Model n300 (with saturation occurring 



only at the highest densities) and nlOOO (which is close to 
fully saturated). The relationship between X and JVh 2 (Fig. 
[6]), however, clearly demonstrates the saturation of the CO 
line. For the models with lower density or metallicity, where 
the W - Nh 2 correlation is rather complex, the X factor can 
have a large distribution, and can even have multiple peaks 
in its PDF (such as Model nlOO, FigJSJi). 




Figure 5. X factor for four models. Nqo is plotted as a function of Nn 2 . The color of each point indicates the X factor. Inset figures 
show the color scale and PDF of the X factor. The corresponding maps of Nn 2 , iVcOi and the X factor from Model n300-Z03 are shown 
in Figure [4] 



4 DISCUSSION AND COMPARISON TO 
OBSERVATIONS 

4.1 The Column Density Distributions: Can CO 
Observations Reveal Log-Normal PDFs? 

Observed morphological characteristics, such as the den- 
sity distribution, may be signatures of the internal dy- 
namics of a cloud. Numerous simulations have shown 
that the (3D) volume density PDF is log-normal for 
non self-gravitating su personic turbulent systems (see 
iMcKee fc Ostrikerll2007l . and references therein) , which may 
be translated into log-normal PDFs of projected (2D ) col- 
umn density (|Ostriker et al.ll200ll ; [Brunt et al.ll2010bl lah . In- 
cluding gas self-gravity, which may counteract the dispersal 
due to turbulence so that some overdense regions may per- 
sist as bound obje cts, results in a PDF with a high density 
tail (|Klessedl200(ih . These signature high density tails have 
been found in column density PDFs obtained fr om observa- 
tions of evolved clouds (|Kainulainen et aUl2009f ). In general, 
however, directly correlating features of the PDF to physical 
processes in the simulation has proven to b e rather difficult 
i|Klesse J2000l ; iFederrath" et al. | |200S| , |2010D . 

Understanding the relationship between CO and the un- 



derlying density distribution is essential if CO observations 
are employed as (molecular) mass tracers of MCs. Usually, 
density distributions are derived from 13 CO intensities, since 
13 CO is optically thin (relative to 12 CO). In our analysis of 
12 CO emission, we find that W does not faithfully trace 
Nco even in some low density regions. This is suggestive 
that 13 CO intensities may also not be an ideal tracer of the 
underlying CO distribution. Though we analyze the density 
distribution traced by 12 CO emission in this work, the suit- 
ability of 13 CO as a column density tracer remains an open 
question that we plan on investigating in future work. 

In practice, the shape of observed density PDFs 
strongly depend s on the tracer observed. For instance, 
iGoodman etail (|2009h found that dust observations of 
the Galactic MC Perseus show log-normal column density 
PDFs, and may provide reasonable estimates of the to- 
tal mass. Gas ( 13 CO) based measures result in PDFs with 
significantly different shapes, though a log-normal fit can 
loosely describe the distribution. In this section, we com- 
pare the PDFs of observed CO intensities with those of the 
intrinsic H2 and CO column density (see Paper I for a dis- 
cussion on the volume density PDFs). 

We have found that the observed W PDF can vary 
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significantly from the column density PDFs, Afco or A/h 2 , 
even though Ah 2 (and Ntot) f° r all models can be well de- 
scribed as log-normal distributions, as shown in Figure [2] 
Only models n300 and nlOOO seemingly have Afco and W 
PDFs that can be well described as log-normals. For the 
PDFs of those models (Fig [2ji,b,e,f), we have overplotted 
best fit log-normals, with functional form 

PDFoce- (loga; - log:co)/2CT2 , (13) 

where x is iVtot, Ah 2 , Afco, or W, a 2 is the variance, and xo 
denotes the mean value. The distributions shown in Figure 
[5] show that the offsets in the log-normals do not coincide. 
Of course, this discrepancy is partly due to our choice of the 
axis ranges of Afco and W. Thus, a comparison of er, which 
is a measure of the width of the log-normal, is the parameter 
best suited for direct comparison between the PDFs. 

For both models n300 and nlOOO, the Atot and Nh 2 
PDFs have a = 0.3. For Model n300, a = 0.5 and 0.2 for the 
Afco and W PDFs, respectively. The best fits give a — 0.3 
and 0.1 for the Afco and W PDFs, respectively, for Model 
nlOOO. As can be seen from the PDFs themselves, the best 
fits Afco and W log-normals are rather different for the two 
models. Such discrepancies are not surprising, given that 



the other models exhibit PDFs of W and Afco with vastly 
different forms. 

That the W and Afco PDFs from Model nlOOO are 
both log-normal, albeit with different characteristics, may 
be somewhat unexpected, given that this model is clearly 
affected by saturation (Figs. [1] and [6}. With line saturation, 
one would expect the W PDF to be asymmetric, with a 
prominent peak towards large intensities, and a sharp drop 
off beyond the peak - the "piled-up" effect discussed in Sec- 
tion 13.11 At any given extinction from Model nlOOO, even 
though much of the line emission from the highest density 
gas is saturated, there is a sufficiently large range of W such 
that there is no sharp drop off in the W PDF beyond the 
peak. This demonstrates that line saturation may be signifi- 
cant even when the PDF of the observed intensities does not 
show a sharp decrease. For such scenarios, it is only when 
inspecting W with AT to t (Fig. [TJ or X with Ah 2 (Fig. 
for which independent estimates of A^tot or Ah 2 would be 
required, that the presence of saturation can be inferred. 

In general, we find that the W PDFs are not log-normal, 
even though the underlying gas column densities are. One 
reason for the difference is that the Afco PDFs themselves 
are not log-normal. Further, in most of the models, the W 
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PDF itself does not neatly follow that of Nco- Since CO 
traces the dense gas, iVco is more intermittent than the un- 
derlying gas distribution. Due to the effects of saturation, 
this intermittency is not captured by W. In the two models 
with log-normal W and iVco PDFs, the variances between 
the distributions are nevertheless different. Taken together, 
we conclude that the distribution of observed CO intensities 
is not an accurate measure of the underlying distribution of 
molecular gas. T his conclusion is affirm ed by the observa- 
tional analysis of lGoodman et all (|2009l ). who found signif- 
icant differences in the 13 CO and dust-based distributions. 



4.2 Globally Averaged X Factor 

We have shown that, in general, the CO integrated intensity 
does not linearly trace Nu 2 , even within individual MCs, 
so that the X factor is not constant. The intrinsic cloud 
properties, such as density, metallicity, and background UV 
radiation field, all influence how W is related to Nu 2 ■ In the 
previous sections, our analysis has dealt with all positions 
in the 2D "observed" maps, which are idealized synthetic 
maps with high spatial resolution of ~ 0.1 pc (the size of the 
grid zone in the MHD simulations). However, in practice the 
derived X factor is likely to be an average over a range in 
extinction values, especially in observations of extragalactic 
clouds where the resolution may be comparable to the cloud 
size. In order for a more direct comparison with these types 
of integrated observations, we now shift our focus to the 
discussion of average quantities in the synthetic maps. 

Figure [7] shows the mean X factor in different Ay bins 
for the four models shown in Figures [5j6j as well as Model 
n300-Z01. The emission-weighted mean X and Ay for each 
model is also shown (by the large symbols). A number of 
features of the averaged X factor are similar to those al- 
ready discussed: a decrease in X with increasing extinction 
for models with low densities or metallicities, the opposite 
trend for Model nlOOO (and to some extent Model n300) due 
to line saturation, and only a small variation in X at high 
densities for Model n300. Qualitatively, the mean values of 
the X factor binned with extinction resemble the extinction 
averaged quantities presented in Paper II (see their Fig. 8). 

One clear difference is that values of the X factor av- 
eraged over a range in extinctions for a given model do not 
span as wide a range as when all positions are considered 
(compare, e.g. to Figs. [5][6]). Only when comparing all the 
models at the low extinction Ay ,$ 5 does the large ex- 
tent in the X factor (10 20 < X < 10 24 cm" 2 K" 1 km" 1 s) 
become clear. 

The striking feature of Figure [7] is the limited variation 
in the X factor at large extinctions, both for a given model 
as well as between models. At extinctions £ 7, X falls in a 
narrow range ~ 10 20 < X < 10 21 cm" 2 KT 1 km" 1 s for all 
models^ This rang e is similar to the X factor values derived 
for the Galaxy JSolomon et al.ll 19871 ; lYoung fe Scovillelll99ll ; 
iDame et al.ll2001 ). The emission weighted mean X and Ay 



values (large symbols) also portray a near constant AGai 
value between the MCs, with the exception of the very low 
metallicity MC. Thus, only considering CO bright regions, 
corresponding to high extinctions, may not expose any real 
variation in the Nn 2 /W ratio, as discussed in Paper II. 

Observations of the Small Magellanic Cloud (SMC), 
which has lower than G alactic metallicity, has resulte d in 
high X factor estimates (|lsraelll 19971 ; iLerov et al.ll2009l . and 
references therein) . Averaged over a whole cloud in the SMC, 
ILerov et al.l |2009) found that the X factor varies by £ 1 
order of magnitude, whereas the X factor computed at the 
CO peaks is similar to the Galactic value. Our analysis of 
CO emission from models with lower metallicity, especially 
Model n300-Z03, which has a metall icity similar to that of 
SMC l| Vermeil fc van der Hulstll2002l ). quantitatively repro- 
duces the estimated X factor, both the larger value due to 
diffuse gas and the lower Galactic value found in the densest 
regions. 



4.3 X Factor in Dense and Diffuse Gas 

In our discussion of the Galactic X factor , we have focused 
on observations of molecu lar clouds (e.g. lYoung fc Scovillel 
Il99ll ; ISolomon et aT1ll987l ) . H owever, CO has also been de- 
tecte d in the diffuse ISM (e.g. iPolk et al.lll988l ; iBurgh et all 
120071 ). and the X factor for these diffuse lines of sight is 
found to be consistent with the stand ard galactic valu e 
of A G ai (few xlO 20 cm -2 K _1 km -1 s. iLiszt et ai]|201Ch . 
Though we consider models of individual (or large volumes 
within) MCs, we investigate whether regions with different 
densities are comparable with the observations of different 
phases of the ISM, focusing o n the recent discus sion on the 
uniformity of the X factor in ILiszt et"af (2010), hereafter 
LPL10. 

LPL10 find that diffuse gas, with molecular column den- 

-2 

D 23 ' 

have X factor values that are ~Xcai- Since the X factor 
can be written as 



X 



(Nh 2 \ = iVco 1 
1 V I \Nco) ~ W /co' 



( Nco \ 
\ W ) 



(14) 



6 This threshold extinction beyond which the X factor is con- 
stant is larger by a factor of ~2 than that found in Paper II. The 
discrepancy arises because the analysis in Paper II relied on sim- 
plified assumptions to compute the CO intensities, whereas here 
we use more accurate radiative transfer calculations. 



LPL10 attribute the constancy of the X factor to the 
increase in the specific brightness, W/Nco, in low den- 
sity, warm gas, by a factor comparable to the decrease in 
/co = Nco /Nn 2 jj For diffus e gas, LPL10 use a n estimate 
of /co « 3 x 10" 6 , based on IBurgh etafl l|2007l ), while for 
dense gas, they adopt /co ~ 10 -4 . These values correspond 
to W/Nco w 10" 15 and 3 x 10~ 17 K km s" 1 cm 2 , in dif- 
fuse and dense gas, respectively, for a constant X factor 
= 3 x 10 20 cm" 2 K" 1 km" 1 s. 

As discussed, we have found that the Milky Way model 
n300 does indeed result in X~Xgh1 in most lines of sight. In 
this model, there is only a very small fraction of the sim- 
ulated volume with the lowest H2 column densities Nn 2 
10 21 cm" 2 (see Fig. [5^0 . Thus, for comparison with 
LPL10 we will use model n300-Z03, which has many lines 
of sight with both dense (Nh 2 ~ 10 22 ) and diffuse gas 



7 LPL10 use the symbol Xqq, whereas we will use /co to avoid 
confusion with the X factor. 



14 R. Shetty et al. 




Figure 7. Mean X factor in Ay bins for 5 models. The X factor is averaged in different Ay bins. The value of X is plotted on the 
midpoint value of Ay of each bin. Each model is identified by different colors and symbols (and labeled in the legend). The large symbols 
shows the global (emission weighted) mean X factor and mean Ay from each model. 



(iV H2 ~ 5 x 10 20 ). Maps of the various relevant quantities 
for this model are presented in Figure [4] Though this model 
has lower metallicity [Z = O.3Z0) than the observations of 
LPL (Z=Zq), a comparison of the observations with Model 
nlOO, (or the high density gas in Model n300) bear similar 
results as those provided by Model n300-Z03. 

We identify dense gas as regions with \og(W) > 
1.5, marked by dashed contours in Fig. [4] and compute 
W/Nco G 4x l(T 17 -8x 1CT 16 K km s" 1 cm 2 . The lower and 
upper limits have corresponding fco values of 5 x 10 -5 and 
5 x 10~ 6 , respectively. The resulting value of the X factor is 
~ 1 - 6 X 10 20 cm" 2 K" 1 km" 1 s. These values of W/N CO , 
fco, and X, reasonably reproduce the values determined by 
LPL10 for dense clouds, with an X factor which only differs 
by at most a factor of 3. 

The low density regions of Model n300-Z03, marked 
roughly by the log(VK) < -1 contour in Fig. [4j have 
W/Nco ~ 10" 15 K km s" 1 cm 2 , and fco ~ 6 x l(r 10 , 
resulting in an X factor ~ 10 cm' " 2 K" 1 km" 1 s. This X 
factor differs from the LPL10 estimated value from obser- 
vations of diffuse gas by up to 4 orders of magnitude. The 
discrepancy can be almost fully attributed to the difference 
m fco between our models and the LPL10 observations; the 
values of W/Nco of our models and that determined by the 
LPL10 analysis only differ by a factor of ~ 2 - 10. Table [2] 
summarizes the comparison of W/Nco, fco, and X between 
Model n300-Z03 and LPL10 in the dense and diffuse gas. 

The comparison of the characteristics of the dense re- 
gions has quantitatively affirmed some level of accuracy of 
the chemical (and MHD) modeling, as well as the radiative 
transfer calculations. Additionally, we can draw some in- 
sights from the discrepancy we find in the characteristics of 
the diffuse regions: the larger difference in fco in our models 
compared to the LPL10 observations of tenuous gas, rather 
than W/Nco, indicates that the chemical processes at work 
in the low density regions within our molecular clouds are 
very different from the chemistry of diffuse Milky Way gas. 
That fco is significantly lower in the low density regions of 



the model suggests that either CO forms at a greater rate in 
the diffuse ISM (relative to the rate in the model), or that 
CO is efficiently transported into the diffuse ISM (and has 
sufficient shielding against photodissociation). The physi- 
cal processes at work in our models, turbulence, magnetic 
fields, or simply even the uniform metallicities (or initial 
densities), may be significantly different in the tenuous ISM. 
In particular, non-thermal chemistry in C-shocks or turbu- 
lent vortices may play a n important role in the chemistry of 
the diffuse gas (see e.g. Federman et all 1 19961 : ISheffer et al.l 
120081 : iGodard et al.l f2009). Alternatively, large scale galac- 
tic processes not considered in our model may be required 
to accurately model the chemistry in diffuse molecular gas, 
such as galactic rotation, the gravitational field due to the 
disk, and/or large scale magnetic fields. We hope to pursue 
studies involving such large scale dynamics and chemistry 
in the future. 



4.4 Future Work: X factors as Mass Estimators 
for CO Bright Cores? 

In observational investigations where the X factor is derived 
solely from CO observations, the CO linewidths are used to 
derive masses through the virial theorem. If the bulk of the 
gas is in molecular form, then W is employed along with 
the CO linewidth derived Nh 2 to estimate the X factor via 
Equation [1] The key assumption h ere is that the observed 
clouds are in virial e quilibrium (e.g. lYoung fc Scoville|[l99ll : 
ISolomon et"aflll987ft . which is the interpretation of the ob- 
serve d power-law l i newidth a v size R rel ationship a v oc R 5 
(e.g. iLarsonl Il98ll : ISolomon et a.1.1 fl987l ') . Measurements of 
the molecular mass through independent observations, such 
as dust-based emission or extinction, can also provide esti- 
mates of Nh 2 ■ These can be combined with CO observations 
to calculate the X factor. As discussed in Section [T] a dis- 
crepancy in the derived X factor based on virialized clouds 
or based on independent H2 mass estimates is found in low 
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Table 2. Comparison of Model n300-Z03 and lLiszt et al.1 (201(1) Observations 



Gas Phase 


AT H2 (cm" 2 ) 


fco 


W/N co 


X (cm" 2 K" 1 km" 1 s) 


Observed Dense Gas (LPL10) 
Dense Gas in Simulation (n300-Z03) 


~ 10 23 
~ 2 - 4 x 10 22 


- 10" 4 

~ 0.5 - 5 x ur s 


~ 3 x 10~" 17 
~ 0.4 - 8 x 10~ 16 


~ 3 x 10 20 
~ 1 - 6 x 10 20 


Observed Diffuse Gas (LPL10) 
Low Density Gas in Simulation (n300-Z03) 


~ 10 20 - 10 21 
~ 2 - 7 x 10 20 


~3x 10~ 6 
~ 6 x 10 -10 


~ io- 15 
~ io- 15 


~3x 10 20 
~ 10 24 



metallicity systems, such as the SMC l|Bolatto et al 



iRubio et al.l(2004llsraej|l997l ; lLerov et al.ll2007l . 2009) 



200S 



The velocity dispersions in our models follow the ob- 
served <r v oc R 5 relationship (Shetty et al. In preparation), 
as found in numerous similar turb u lent cloud models (e. | 
iKlessenl |200d : lOstriker et all l200ll ; iFederrath et~ai1 I201C 
We have calculated the X factors in model MCs assum- 
ing independent knowledge of Ah 2 , s ° no assumption of a 
linewidth-size scaling, or of virialized clouds (or cores), is 
necessary. In our analysis of total column densities and in- 
tegrated intensities, we find that X factor is primarily con- 
trolled by the CO abundance, and thus ultimately by cloud 
properties most responsible for CO formation: the metallic- 
ity, density, or background UV radiation field. This variation 
in the X factor is most drastic - upto several orders of mag- 
nitude - in models with metallicities, densities, or UV radi- 
ation fields different from the Milky Way. These results are 
in agreement with models s howing an increased X factor 
for lo w metallicity systems l|Malonev fe Blackl 1 19881 : llsraell 
1 19971) . 

The discrepancy between X factor estimates based on 
virialized systems, and those based on independent Ah 2 
measures can be explained by the selective photodissociation 
of CO in low density regions. The CO bright objects thus 
only trace highest density gas, which is surrounded by lower 
density molecular material, with little or no CO emission 
ilMalonev fc Blackl I l988l; iBolatto et al.ll2008l : iGrenier et all 



2005!: IWolfire et al 



2010, Paper II, Molina et al. in prep. 



In low metallicity or low density systems, we only find the X 
factor to be constant at the very highest densities (Figs. |6p 
&[6ji), though the vast majority of observed points clearly 
show a gradient in the X factor with density. For the high 
density Model nlOOO (Fig[6j;), the X factor increases with in- 
creasing density due to CO line saturation (Fig. []J) , though 
the range of X is only £ 2 orders of magnitude. 

The idea that CO observations accurately provide mass 
estimates, with XmXaax, of virialized clouds needs to be 
affirmed quantitatively. Previous efforts t o address this is- 
sue h ave usually involved static models. iKutner fc Leunel 
i|l985l) constructed clouds models with microturbulent ve- 
locities consistent with the observed linewidth-size relation- 
ship. They found that the X factor still strongly depends on 
the temperature, as well as the density and CO abundance . 
Using photodissociation region models, IWolfire et ail (|l993l ) 
found that microturbulence cannot reproduce the observed 
CO line profiles. The CO intensity, and hence the X factor, is 
sensitive to other parameters responsible for CO formation. 
With the radiative transfer calculations performed on the 
3D MHD simulations including chemistry discussed here, we 
can further investigate how various turbulent velocity fields 



influence CO emission, as well as the cloud characteristics 
that affect the X factor in a range of environments. 

Such an analysis would address the following is- 
sues: How well do CO peaks trace coherent objects 
in spectral (PPV) cubes? Identifying coherent objects 
from spectral cubes has proved to be quite chal- 
lenging, as a conseq u ence of line-of-sight projection 
i|Adler fc Roberts! 1 19921 : iPichardo etHI |2000D . How well 
do CO linewidths encode information on the dynamics of 
the cloud? Projection effec ts may also skew the linewidth 
- size power law scaling s ( Ballcs teros-Paredes et all 1 19991 : 
iBallesteros-Paredes fc Mac Lowll2002l ; IShettv et al.ll201oC lf 
a CO bright core can be accurately identified, is the observed 
linewidth representative of the intrinsic velocity dispersion? 
To what extent does turbulence affect the CO emission, and 
thus the X factor? These issues will all be visited and dis- 
cussed in a forthcoming paper. 



5 SUMMARY 

We have performed radiative transfer calculations to inves- 
tigate the nature of 12 CO (J=l-0) emission from simulations 
of molecular clouds (MCs). The MCs are modeled through 
hydrodynamic simulations of a turbulent, magnetized, non 
self-gravitating gaseous medium, along with a treatment of 
chemistry to track the formation of H2 and CO (Paper I 
and Paper II). As part of the radiative transfer calculations, 
we use the Sobolev (LVG) method to solve for the CO level 
populations. We analyze the probability distribution func- 
tions, and X factor properties, using the velocity integrated 
CO intensities W , along with the column density of CO and 
H2, Aco and Nn 2 , respectively. Our main findings are: 

1) In all models, W increases with increasing total hy- 
drogen column density N to t (or extinction Av). However, for 
the Milky Way and high density models (n300 and nlOOO), 
which have the highest CO abundances, we find a threshold 
in W«65 K km s~ at high extinction due to saturation of 
the CO line. 

2) All models have log-normal ATtot and Nn 2 distribu- 
tions. In general, however, the W and Aco PDFs are not log- 
normal. Further, since CO is optically thick, the W PDFs 
do not have similar shapes to the corresponding Aco PDFs. 

3) In some models for which the CO line is saturated, 
the peak in Aco is offset towards higher densities than the 
peak in W, though the two PDFs seem to be correlated at 
low densities. However, such a "piled-up" W PDF does not 
necessarily arise in clouds with saturated CO emission, espe- 
cially those with a limited range in densities (as for the high 
density model nlOOO, Figs. [T][2]). Independent measurements 
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of N to t or Nh 2 , along with W , are needed to unambiguously 
identify saturation. 

4) The X factor is not constant within individual molec- 
ular clouds, and in models with low CO fractions, can vary 
by up to 4 orders of magnitude. The low density regions have 
the highest X factor, in agreement with previous modeling 
efforts. 

5) In most simulations, the averaged X factor is found 
to be similar to XqsI (~ 2 x 10 20 cm -2 K _1 km -1 s), mask- 
ing the variation of X with density within clouds. In clouds 
with low CO abundances relative to the Galaxy (such as the 
model MC in a dwarf galaxy or the LMC/SMC), the densest 
regions have X ~ -Xciai- 

6) Emission weighted averaged X factors from all mod- 
els provide values ~XGai with the exception of the very low 
metallicity MC (n300-Z01). Similarly, at extinctions Av ^ 
7, the X factor for always falls in a narrow range ~ 10 20 
£ X < 10 21 cm -2 K" 1 km -1 s. As discussed in Paper II, 
observations targeting CO bright regions may be unable to 
detect real variations in the X factor. 

7) In the low density gas within the MC models, the 
CO fraction fco is found to be ~ three orders of ma gni- 
tude lower than fco measured by iBurgh et al l |2007h in 
diffuse gas. That -/Vh 2 and W /Nco are comparable between 
the model and the observations suggest that the dynamics 
and chemical evolution of diffuse gas in our molecular cloud 
model are rather different from that in the large scale diffuse 
gas of the Galaxy. 

8) We do not assume that clouds or cores are in virial 
equilibrium, and the X factor variations we find are in gen- 
eral agreement with observational analyses employing inde- 
pendent measures of total molecular column densities. In a 
follow up investigation, we will use the spectral information 
of our models to assess the X factor when virial equilibrium 
is assumed; we will then be able to address the X factor 
discrepancy found in observational works when assuming 
virialized CO clouds or when using independent molecular 
mass estimates. 
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